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The viscosity of aeolotropic liquids between two plates has been found to increase when an 
electric field, steady or alternating, is applied perpendicularly to the plates; this effect is probably 
due to real or apparent orientation of the axis of the swarms parallel to the electric field. A 
change of electrical conductivity under shearing stress has also been found. Ordinary isotropic 
insulating liquids show some apparent decrease in viscosity with a steady field, but this is 


probably due to electrolysis. 


HE liquids called mesomorphic melts, liquid 
crystals or aeolotropic liquids display a 
number of anomalies in comparison with ordinary 
liquids. Of the most striking of these anomalies 
we mention the following. , 
Viscosity shows a peculiar variation with the 
temperature. The phenomenon was studied at an 
early stage in the history of the mesomorphic 
liquids and was used in order to disprove the 
emulsion theory. ; 
The electrical conductivity due to the presence 
of dissolved electrolyte is changed by the action 
of a magnetic field, a phenomenon discovered by 
The. Svedberg.! The influence of a magnetic or an 
electric field on the light. absorption has been 
studied by the author.? Even mechanical action 
changes the light absorption.? The influence of a 
magnetic field on the diffusion constant of a 
solute has been investigated by Svedberg.’ These 
effects are only found in liquids of the nematic 
type. For several other phenomena I refer to the 
summary published by Oseen.* 

!The. Svedberg, Uber die Elektrizitatsleitung in aniso- 

tropen Flissigkeiten, Ann. d. Physik. 44, 1121 (1914). 
* Y. Bjérnstahl, Untersuchungen iiber anisotrope Flissig- 


keiten, Ann. d. Physik 56, 161 (1918). 

* The. Svedberg, Diffusion in anisotropen Fliissigkeiten. 
Kolloid Zeits. 22, 68 (1918). 

*C. W. Oseen, Die anisotropen Fliissigkeiten (Fortschritte 
Physik und physikalischen Chemie, 20, 

0. 2). 


It is therefore perhaps not astonishing that the 
viscosity of certain liquid crystals of the nematic 
type shows a marked dependence on an electric 
field. As regards ordinary liquids it has generally 
been assumed that such an action does not exist. 
This view is in accordance with theoretical con- 
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Fic. 1. Rotor. N, clamp; M, mirror; Q, disk; O, weight. 
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Fic. 2. Schematic diagram of the apparatus. R, rotor; 
P, electrode; M, mirror; W, heater; L, stirrer; B, case; H, 
torsion head; G, pointer; C, chuck; D, micrometer head; 
A, insulator; K, telescope; J, J insulators. 


siderations.° I will show, however, that there are 
certain apparent exceptions to this rule. 


EXPERIMENTAL ARRANGEMENT 


The apparatus consisted in principle of a rotor, 
a disk suspended by a wire and immersed in the 
liquid to be investigated, and a fixed plate below 
the disk. The two plates formed the two termi- 
nals. The damped oscillations of the rotor were 
studied. 

The rotor was made of brass (Fig. 1). A clamp 
N and a mirror M were fixed at the top. The 
weight O was added in order to increase the 
moment of inertia. The disk Q was gold-plated 
in some cases. 

A sectional drawing of the whole arrangement 
is given in Fig. 2. The rotor R was suspended by a 


Gerlach, Materie, Elektrizitaét, Energie, Dresden 
(1926), p. 40. 
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phosphor bronze wire, fixed in the chuck C, which 
was attached to the torsion head /7. One part of 
this could be rotated and was provided with a 
pointer G. The screw E could be moved up and 
down by means of the nut F. The position of the 
pointer G could be read on a dial #7. The torsion 
head was insulated from the rod of the support by 
means of the insulator A and from the case by 
means of J. The distance between the rotor disk 
and the plate P was determined with the aid of 
the micrometer head D, on which 0.01 mm could 
be read. The reading was checked by means of a 
telescope K, which was provided with a micro- 
metric eyepiece. 

The liquid was contained in a cylindrical vessel 
of Jena Normal Glass, at the bottom of which a 
brass plate P was fixed by means of a platinum 
bolt. This was connected with the binding post 
P». The distance between the upper surface of 
the plate P and the free surface of the liquid was 
20 mm. The glass cylinder was immersed in paraf- 
fin oil, contained in a Dewar beaker. The oil was 
heated electrically by means of a resistance W 
and efficient circulation was obtained with the 
stirrer L. The temperature of the oil was deter- 
mined with the aid of a thermometer graduated 
into 0.2°C. The difference between the tempera- 
ture of the liquid in the tube and that of the oil 
was measured with a thermoelement. 

The amplitudes of the rotor were measured 
with the aid of mirror, scale and telescope. In 
order to avoid corrections for large deviations, a 
circular scale of 648 mm diameter was used in 
the center of which the rotor was placed. The case 
(B) prevented air currents from disturbing the 
movements of the rotor. A glass window at the 
bottom of the case permitted observation of the 
oscillations of the rotor. The time of oscillation in 
air was 25.8 sec. 


Electrical connections 


A source of direct current was connected to the 
terminals of a potentiometer, from which the 
required voltage was obtained. 

Alternating current was obtained from the 
secondary of a transformer. The terminals of the 
primary were connected to an inductor alterna- 
tor, the voltage of which could be controlled. The 
frequency was 1200 cycles per sec. The voltage 
was measured by means of a hot-wire voltmeter. 
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Device for filling and emptying the tube 


The substances investigated, e.g., p-azoxy- 
phenetol, are solid at room temperature. p-azoxy- 
phenetol melts and passes into the mesomorphic 
state at a temperature of 137°C. The liquid be- 
comes isotropic at 167.5°C. Even the meso- 
morphic liquid has a considerable vapor tension. 
In order to fill and empty the tube and to clean 
it easily without having to change its position, 
operations which are often carried out, a special 
device was designed. 


The size of the rotor and the disk 


It is not easy to find a suitable size for the 
rotor. The time of oscillation should not be less 
than about 10 seconds, since it is otherwise diffi- 
cult to take the readings. The time of oscillation 
can be increased by increasing the moment of 
inertia by a weight. The diameter of this weight 
is limited by that of the tube, which must be 
cylindrical. The weight of the whole rotor must 
not exceed a value corresponding to the elastic 
limit of the wire. Finally the logarithmic decre- 
ment must possess an adequate value. 

Several different types of rotors were designed 
and constructed. In order to check the constancy 
of the torsion modulus, the time of oscillation 
was determined after every experiment. 


VISCOSITY OF AN AEOLOTROPIC LIQUID 
DEPENDS ON THE MAGNITUDE OF THE 
SHEARING STRESS 


A comparison of the values obtained by using 
our apparatus with those obtained by the capil- 
lary method may be of interest. 

The relation between the logarithmic decre- 
ment A, the viscosity.» and the density p of a 
liquid can be expressed by an equation, derived 
by Meyer® and used by Fawsitt? and others for 
the determination of the viscosity of melts. The 
equation is, 


A—Ao= Ci(pn)'+ Copn+ C3n, 


where Xo is the logarithmic decrement in air, and 


my ‘in 6). Meyer, Innere Reibung, Poggendorfs Ann. 113, 
7C. H. Fawsitt, On the Determination of Viscosity at high 
Temperatures, Proc. Roy. Soc. A80, 290 (1908). 
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Fic. 3. Influence of the experimental arrangement on the 
viscosity. Closed circles, values obtained by Eichwald and 
published by Schenck.® ‘Open circles, values obtained with 
the aid of the rotor method. 


C;, C2 and C; are constants. As will be shown be- 
low this equation is not strictly applicable to the 
experimental arrangement in question. An in- 
vestigation has shown that the equation, after 
determination of the constants C; can be used as 
an interpolation formula. 

The results of some measurements on p-azo- 
xyanisol are recorded in Fig. 3. The values on the 
viscosity of the isotropic liquid are identical with 
those obtained by the capillary method. This, 
however, is not the case for the liquid in the 
aeolotropic state. The values calculated from 
experiments with the rotor method are greater 
than those obtained by Ostwald’s method. The 
reason for this is that a more complete orienta- 
tion of the particles or a more complete deforma- 
tion of the liquid is effected when the last method 
is used. The result is of interest since it shows that 
the value of the viscosity is a function of the 
shearing stress which in turn is dependent on the 
experimental arrangement used and on the 
structure formed in the liquid. 


_=R. Schenck, Kristallinische Flissigkeiten, (1905) Leip- 
zig. 
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THE INFLUENCE OF AN ELECTRIC FIELD ON 
Isotropic LIQUIDS 


With the aid of the arrangement described 
some preliminary experiments were made in order 
to ascertain if ordinary isotropic liquids are in- 
fluenced by an electric field. Table I shows some 
results with xylol and direct voltage. The dis- 
tance between the fixed plate and the disk was 
2.0 mm. In the table, V the voltage is in volts, T 
the temperature and \ the decrement. 

If the rotor is the positive terminal a decrease 
of decrement is found, if it is negative an increase 
results. The latter effect is more marked. The 
experiment was repeated with alternating volt- 
age. Some values are given in Table II. 


TABLE I. Influence of an electric _ xylol. (Time of 


oscillation for V =0, 26.0 sec.) 
V T 
0 19.2° 0.0270 
+220 re 0.0261 
—220 0.0290 


TABLE II. Influence of an electric field on xylol. (Time of 
oscillation 26.0 sec., V=0). 


V T 
0 20.2 0.0261 
150 0.0264 
0 0.0267 
150 ” 0.0264 


TABLE III. Influence of an electric field on isotropic 
a Distance between plates, 2.50 mm. Time 
of oscillation for V=0, 26.2 seconds. The rotor was the 
negative terminal. JT is the temperature in °C, V the 
voltage in volts. 


T V/cm 
179.1 0 0.068 
180.0 160 0.062 
180.0 320 0.053 
179.8 0 0.067 
179.7 640 0.045 


TABLE IV. The change of decrement by an electric field. 
Isotropic p-azoxyphenetol, alternating voltage. 


V/cm 
174.0 0 0.074 
174.0 400 0.075 
174.0 0 0.072 
174.1 600 0.073 
174.2 0 0.075 
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Fic. 4. Dependence of the field on the logarithmic decrement 
(isotropic liquid; direct voltage). 


In this case no change in decrement can be 
found. Similar results were obtained with other 
liquids, for instance with benzene. Some values 
from experiments on azoxyphenetol in the iso- 
tropic state are recorded in Table III, and dia- 
gram 4. 

The decrement is decreased with the voltage. 
The decrease was less when the rotor was the 
positive terminal. 

The results of some experiments with alternat- 
ing voltage are presented in Table IV. 

These measurements show that the decrement 
is independent of the voltage of the alternating 
current. 

The effect found with direct voltage can 
probably be interpreted as due to an electrolytic 
action. p-azoxyphenetol possesses a conductivity 
of the order of magnitude 10-* ohm~'cm~. The 
surface layer of the oscillating disk was changed. 
A change of the viscosity of the liquid, ie., a 
volume effect seems to be excluded at least in 
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Fic. 5. Dependence of an electric field on the decrement 
(direct voltage). 
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° alfernating voltage 
@ direct voltage 
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0.03 
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Fic. 6. Dependence of temperature on the decrement. 


most cases. The conclusions regarding the orien- 
tation of the molecules that can be drawn from 
Langevin’s theory for the Kerr effect show that 
even in most liquids with a marked dipole mo- 
ment this orientation is insignificant. 


THE INFLUENCE OF THE ELECTRIC FIELD ON THE 
AEOLOTROPIC LIQUID 


Some results obtained from measurements of 
the change of the logarithmic decrement in an 


TABLE V. Change of decrement in an electric field. 


(p-azoxyphenetol) direct voltage, time of oscillation 
26.2 sec. (V=0). 


T V/em p 
165.0 0 0.065 1.070 0 
165.0 40 0.067 0.002 
165.0 160 0.079 0.014 
165.6 80 0.072 0.007 
164.7 120 0.074 0.009 
164.3 240 0.085 0.020 
164.7 480 0.091 0.026 
164.7 800 0.088 0.023 
155.0 0 0.075 1.082 0 
154.2 160 0.089 0.014 
154.3 320 0.111 0.036 
154.6 80 0.077 . 0.002 
154.6 120 0.087 0.012 
154.5 240 0.104 0.029 
154.7 480 0.114 0.039 
155.0 800 0.114 0.039 
145.1 40 0.084 1.092 
146.2 104 0.094 
145.4 240 0.119 
145.3 480 0.129 
145.4 800 0.127 
141.1 0 0.083 1.096 0 
141.4 40 0.081 —0.002 
141.4 80 0.088 0.005 
141.4 120 0.098 0.015 
141.4 240 0.128 0.045 
141.4 480 0.146 0.063 
141.4 800 0.140 0.057 
141.6 0 0.080 
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Fic. 7. Dependence of an electric field on the decrement 
(alternating voltage). 


electric field (direct voltage) are given in 
Table V. In Table V, T the temperature is in °C, 
V/cm the intensity of the field in volts per cm, 
\ the decrement and p the specific gravity of the 
liquid. 

TABLE VI. Influence of an electric field (alternating 


voltage) on decrement in p-azoxyphenctol. Distance between 
plates 2.5 mm; time of oscillation 26.0 sec. (V =0). 


T Viem p 
143.7 0 0.078 1.094 0 
143.2 0 0.079 0.00 
143.2 80 0.079 0.001 
143.2 120 0.083 0.005 
143.0 160 0.095 0.017 
143.0 200 0.101 0.023 
143.5 0 0.077 —0.001 
143.5 240 0.113 0.035 
143.3 320 0.121 0.043 
143.3 0 0.079 0 
143.3 400 0.125 0.046 
143.3 480 0.112 0.033 
143.3 600 0.098 0.020 
143.3 0 0.079 0.000 
154.4 0 0.068 1.082 0 
154.9 40 0.067 —0.001 
154.9 80 0.071 0.003 
154.9 160 0.077 0.009 
154.6 0 0.069 0.001 
154.6 240 0.096 0.028 
154.6 320 0.100 0.032 
154.5 0 0.068 0.C00 
154.5 480 0.094 0.026 
154.5 600 0.083 0.015 
154.7 0 0.068 0.000 
165.5 0 1.070 
165.5 40 0.064 0.000 
165.0 80 0.065 0.001 
165.2 160 0974 0.010 
165.2 0 0.064 0 
165.2 240 0.080 0.016 
165.4 320 0.086 0.022 
165.4 480 0.087 0.023 
165.5 600 0.082 0.018 
165.5 0 0.064 0.000 
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The rotor was the negative terminal. The effect 
was independent of the direction of the field. As 
Fig. 5 shows, the decrement is independent of the 
voltage up to about 100 volts per cm. At higher 
voltages the decrement rises and reaches a value 
of saturation at about 500 volts per cm. The 
action is most pronounced at low temperatures. 
The temperature variation of the effect is shown 
in Fig. 6. 

Similar experiments were made with alternat- 
ing voltage. The results are given in Table VI. 

The symbols are the same as in Table V. Fig. 7 
shows the change of the decrement with the in- 
tensity of the field. In Fig. 6 the influence of the 
temperature is shown. The change increases with 
the field, attains a maximum and decreases again. 
The effect is greatest at low temperatures. 

Similar results were obtained with p-azoxyani- 
sol, although the effect was smaller. 


THEORY 


The problem of a disk oscillating in a viscous 
medium has been treated by Coulomb,’ Maxwell" 
and Meyer.’ Coulomb and Meyer discuss the 
case of a disk oscillating in a very large volume, 
while Maxwell assumes that the disk is placed in 
the middle of two fixed plates. For the most part 
the method has been used by these authors to 
study the viscosity of air. 

A numerical lower limit for the change of vis- 
cosity in an electric field can be obtained if we 
assume that the change is produced in the whole 
liquid. A value of ». is obtained from equation 
(p. 5) developed from Meyer’s theory. The values 
obtained using such a theory are of course smaller 
than the actual ones. 

Neither theory can be applied to the experi- 
mental arrangement used by us. However, it is 
possible to generalize Maxwell’s discussion 
and obtain a relation between decrement and 
viscosity. 

Consider a layer of liquid of the height dy at 
the distance y from the fixed plate (Fig. 8) and 
an element of volume at the distance r from the 
axis and with an azimuth ®. If the viscosity of the 
liquid between the disk and the plate be 7, and 


®Ch. Coulomb, Mem. de I’Inst. 4, 431 (1822). 
10 J. C. Maxwell, Proc. Roy. Soc. A156, 259 (1866). 
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the tangential velocity of the element », the force 
which is acting on the lower surface of the volume 
element is ndv/dy or nrd°0/dtdy. The force acting 
on the upper surface of the layer is 


nrd?0/dtdy+ nr(d°0/dy?dt)dy. 


The resultant force acting on the element is 
nrd°0/dy’dt and the equation of motion of the 
element is expressed by 


pd70/dt? = ; (B) 


For a layer between the disk and the surface of 
the liquid at the distance y’ from the surface at 
the liquid we get 


926’ /dt? = /ay’2at, (A) 


where 7’ is the viscosity and 6’ the azimuth of the 
volume element. I have assumed two different 
viscosities » and n’. As a matter of fact the vis- 
cosity in the aeolotropic liquid is not a constant 
but varies during an oscillation and with y and y’, 
respectively, since it depends on dv/dy. Even in 
the same layer the viscosity decreases with in- 
creasing 7, and it is only in the interest of sim- 
plicity that I have made a provisional assumption 
of only two different viscosities 7 and 7’. 

The equations A and B show that the motion is 
independent of r. The boundary conditions are: 


0°70’ /dy’dt=0 fy=0 
™™ 0=0'. 
- The equation of movement of the rotor is ex- 
pressed by 
076 dé 070 076’ 
—+2h—+070 =0, 
at? ot dydt dy’ dt 
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where I is the movement of inertia, A=r*/2 
and the second term is due to the damping action 


of the air on the upper part of the rotor. Let us 
try the solution 


= elie) t+a’ cos u'y’ = eatta’ cos u'y! 
sin uy; 


where u, u’ are complex quantities. From (A) 
and (B) we obtain 


2 
pa=—7'u" pa=—nu?; 


The conditions for y=b; and y’=—a give 
cos u’a= sin ub; 


+ nau ctg ub—n’ au’ tg u’a=0, 


I 2h 1 1 

+—+— —-ctg ub+— tg w'a=0. 
pA a as u u' 

The quantities J, A, h, Q, p, b, a are all known. 

The equation may be written 


n, ]=0 


or if the real and imaginary parts are separated 
and if a= —/+iw where / is the damping factor 
and w the frequency 


¢-Ll, ®, 7, n’ 
log!’ e. 


The quantities / and w are determined by experi- 
ment. We can obtain values on 7 and 7’. It can 
be shown that a numerical calculation can be 
carried out at least if the difference 7—7’ is not 
great.* 

In the case of an electrical field we obtain the 
analogous equations 


drLley Wey Ney ne 
Wes Ney ne 


where 7, is the viscosity under the influence of the 
field, 7.’ the viscosity above the disk, /, the 
damping factor, and w, the frequency at the 


oscillation. Since n,’=7' it is not necessary to 
determine w,. 


*The author is indebted to Dr. E. Hogner and to Mr. 


N. Larsson for some suggestions as regards these calcu- 
lations. 
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Although it is possible to calculate n, n’ and n,. 
from the equations, we refrain from giving any 
numerical values. The theoretical assumptions are 
doubtful. The problem is much more complicated 
than we have assumed. Furthermore, the aeolo- 
tropic liquids possess a special structure, which is 
distorted by the oscillations of the disk. As 
another argument may be mentioned that New- 
ton’s assumptions for the tensor of viscosity is not 
valid for aeolotropic liquids. Oseen™ and Anze- 
lius” have suggested still other possibilities. 


Change of electrical conductivity under shearing 
stress 


As regards the cause of this electrical viscosity 
effect it is not difficult to explain it as resulting 
from an orientation or a distortion of the swarms. 
This hypothesis has been verified by comparison 
with the results obtained from measurements of 
the electrical conductivity and dielectric constant. 

The commutator method used by The Sved- 
berg! for determinations of the conductivity of 
comparatively concentrated electrolyte solutions 
fails in the case of pure liquids, which are char- 
acterized by a small conductivity. However, if 
special precautions are observed, the ordinary 
bridge method can be adapted even to the meas- 
urement of the high resistances in question. An 
earthing device analogous to that designed by 
Wagner" was used and certain of the resistances 
were shielded. Voltage of a frequency of about 
1000 cycles per sec. was obtained from an alter- 
nator, which allowed the voltage to be varied be- 
tween the desired limits. A band-pass filter served 
to suppress the harmonies. In order to prevent a 
too assymetrical distribution of the earth-admit- 
tances an input transformer was used. As null 
indicator an impedance matching transformer 
connected to a two-stage amplifier and a tele- 
phone was used.“ 

Since disturbances due to the capacity to earth 
are possible, the arrangement was tested by the 
following method. A resistance was substituted 
for the cell and the voltage was varied. The re- 


1 C, W. Oseen, Trans. Faraday Soc. 29, 895 (1933). 

2A. Anzelius, Uber die Bewegung der anisotropen 
Flissigkeiten, Diss. Upsala 1931. 

1K. W. Wagner, Elektrotech. Zeits. 32, 1001 (1931). 


4 The author is indebted to Professor L. Ramberg for 
the loan of some instruments. 
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TABLE VII. Dependence of conductivity on the field, 
(p-azoxyphenetol; T = 144°). 


E 108 
(volts per cm) 
6.8 4.0 
19 4.0 
65 4.1 
157 4.6 
195 4.7 


sistance was found to be independent of the 
voltage. 

Since the method was to be used only for the 
purpose of analysis the cell constant was not 
determined, but the specific conductivity was 
calculated from the diameter and distance of the 
plates neglecting boundary corrections. The de- 
pendence of the conductivity on the field was 
determined (Table VII). The conductivity in- 
creases with the intensity. It is evident that the 
critical limit is about 70 volts per cm. In order to 
eliminate this influence a low voltage was used. 

When the rotor was at rest a conductivity 
x=4.5X10-* was obtained. With 
oscillation, i.e., under shearing stress x was found 
to be 4.2 10-8 at a distance between the plates 
1.02 mm, the voltage being 0.67 volts, correspond- 
ing to 6.6 volts per cm. The temperature was 
144°C. In other experiments the decrease of con- 
ductivity varied between 0.0—0.210-* 
em", 

The measurements have been of a preliminary 
nature and are not yet completed. A comparison 
with the results obtained by Svedberg! shows 
that in the absence of an electric field the swarms 
in the liquid under shearing stress display a 
tendency to arrange themselves with their axes 
parallel to the plates. In an electric field an 
orientation is effected with the axis at right angles 
to the plates. 

It is evident that layers consisting of swarms 
orientated in this way must tend to increase the 
viscosity. 
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It has been supposed that the axis of the 
swarms in a nematic liquid tends towards a 
transversal position in an electric field for the 
following reason. In a magnetic field a decrease of 
the dielectric constant has been observed sug- 
gesting that e,<e., where e¢, is the dielectric 
constant parallel to the axis of the swarm. 

In an electric field the stable equilibrium cor- 
responding to a minimum of potential energy 
would thus involve a transversal orientation of 
the axes of the swarms. 

My results militate against the theory of trans- 
versal orientation, but agree closely with those 
obtained by Kast," who has found by x-ray in- 
vestigations that the axis of the swarm inclines 
towards a longitudinal position in an electric 
field. In order to explain this curious phenomenon 
Ornstein'® has suggested that the swarms have a 
fixed dipole moment in the direction of the axis, 
which only slightly affects the value of the di- 
electric constant. 

In view of the chemical formula for p-azoxy- 
phenetol, it seems unlikely that a great dipole 
moment should exist in a direction parallel with 
the axis of the molecule. I refer to my alternative 
hypothesis, first evolved in 1918 concerning the 
orientation of the swarms in a nematic liquid? 
p. 183 and p. 201. 

By microscopic observations on p-azoxyphene- 
tol and p-azoxyanisol in an homogeneous electric 
field, I have found small spherical particles or 
drops, which move about in the liquid. The phe- 
nomenon is probably due to the formation of 
space charges. The shearing stress produced by the 
movement causes the orientation of the swarms. 
That such ‘“‘mechanical”’ orientation is possible is 
clearly shown by the experiment described 
above. The movement of the rotor causes the 
axis of swarms to become parallel with the plates. 


1 W. Kast, Anitsotrope Fliissigkeiten im elektrischen 
Felde, Zeits. f. Physik 71, 39 (1931). 
1 L. S. Ornstein, Zeits. f. Krist. 79, 331 (1931). 
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The Motion of a Viscous Fluid Under a Surface Load 


N. A. HASKELL, Massachusetts Institute of Technology 
(Received April 25, 1935) 


A formal solution is given for the motion of a highly viscous fluid when a symmetrical pressure 
is applied at the surface. This is applied to the subsidence of a cylindrical body of constant 
thickness and to the recovery of the fluid after removal of a load. Applying the latter case to 
the plastic recoil of the earth after the disappearance of the Pleistocene ice sheets, it is found 
that the geological data imply a kinematic viscosity of the order of 3 X10”! c.g.s. units. 


INTRODUCTION 


BSERVATIONS on the elevated and tilted 
shore lines of the glaciated tracts of Europe 
and North America have shown that during and 
following the disappearance of the Pleistocene 
ice sheets the earth returned to a configuration of 
equilibrium by delayed plastic flow as well as by 
immediate elastic recovery. The progress of this 
recovery has been dated by means of varved 
glacial clays and F. Nansen! has constructed 
curves illustrating the uplift of the Fennoscan- 
dian region from 16,000 B.C. to the present time. 
This suggests the possibility of determining the 
effective viscosity of the earth’s outer shells under 
forces of long duration. 

In the present paper we shall neglect the curva- 
ture of the earth and treat the problem of the 
motion of a semi-infinite, incompressible, viscous 
fluid under the action of a radially symmetrical 
pressure applied at the free surface. Since, in the 
case of the earth, we are dealing with extremely 
small accelerations and very high viscosity, we 
may neglect the inertial terms in the equations of 
motion in comparison with those arising from 
viscous forces. We shall first treat the case of an 
arbitrary symmetrical impressed pressure and 
then apply the results to two special cases, (1) a 
constant load of radius ro applied at ‘=0, and 
(2) recovery after the removal of an arbitrary 
symmetrical load. 


1. FORMULATION OF THE PROBLEM 
Neglecting the terms arising from the accelera- 


tion, the equations of motion of a fluid of vis- 
cosity 7 and density p in a gravitational field g 
are ;— 
nv? V=grad p—(0, 0, pg), (1.1) 
div V=0, (1.2) 


1F. Nansen, The Earth's Crust, Its Surface Forms and 
Isostatic Adjustment (1928), Oslo. 


where the positive z axis has been taken as 
directed downward. Transforming to cylindrical _ 
coordinates (r, z, ¢), assuming radial symmetry, 
and setting p= p—pgz, these become 


(1.11) 
(1.12) 
(1.21) 


The components of stress in which we are inter- 
ested are 


= (p+ pgz) +2n0d V,/d2, (1.31) 
(1.32) 


The boundary conditions are that on the free 
surface p,, shall vanish and p,, shall equal the 
applied stress, and at infinity the stresses and 
velocities shall vanish. 

Let the equation of the free surface be z= ¢(r, #) 
and take the undisturbed surface as the plane 
z=0. If we assume that ¢ remains small in com- 
parison with other distances entering into the 
problem, such as the radius of the applied load, 
we may replace the value of dV,/dz at z= ¢ by its 
value at z=0, and similarly with all the other 
quantities appearing in (1.31) and (1.32) except 
pgz.? Calling the applied pressure —o(r, the 
boundary conditions become 


D(r, 0, t)+ pgs (r, t) 


—2n(dV./dz)(r, 0, t)=o(r, t), (1.41) 
(V,/d2+0V,/dr) =0 (1.42) 
and 0¢/dt= V,(r, 0, (1.43) 


2. GENERAL SOLUTION FOR ARBITRARY Im- 
PRESSED LOAD 

Setting V,=Ri(r)Z1(z), Vz=R2(r)Z2(z), p 

=R;(r)Z3(z) in (1.11), (1.12), and (1.21), and 

? This is the sort of approximation commonly made in 


treating waves on the surface of a fluid. Cf., H. Lamb, 
Hydrodynamics, fifth edition, Chaps. IX. & XI. 
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separating variables, we have the Bessel equa- 
tions, 


d/rdr(rdR;/dr) —R,/r?+?Ri=0, (2.11) 
d/rdr(rdR2/dr)+?R2=0 (2.12) 


and R;=const. X Re. (2.13) 


The appropriate solutions are J;(Ar), Rs 
=J,(Ar), where we have included the constant 
factor in R; in the z factor. The z equations are 


—Z3/n, (2.21) 
—\*Z2=dZ;/ndz, (2.22) 
+dZ2/dz=0. (2.23) 


Eliminating Z; and Z;, we have the fourth order 
equation for Z2 


(2.3) 


The solutions of this equation are e*™, se** of 
which only those with negative exponents are 
appropriate to the present problem. The solutions 
of the set (2.21), (2.22), (2.23) are then, 


Z2=e-*(A+Bz), (2.41) 
Z,=e-*(A —B/A+Bz), (2.42) 
Z3=2nBe-. (2.43) 


In order to satisfy the boundary condition 
(1.42) we have B=)A.* We must now satisfy 
(1.41) with functions of the form 


(2.51) 
Vi= fo A (2.52) 
p=2nfo°A (2.53) 


From (2.51) dV,/dz=0 at z=0, hence (1.41) 
becomes 


2nfo°A (dA) Jo(Ar)AdA+ pgf=a(r, t). (2.61) 


A will evidently have to be a function of the time 
in order to satisfy this equation, hence we may 
differentiate with respect to ¢ and use (1.43). 
Then 


So? Jo(ar) {2ndA /dt+ (2.62) 


3 Note that this is the same relationship that we would 
get if instead of the vanishing of p,, at z=0 we took asa 
boundary condition that V,=0 at z=0. The latter is the 
condition we would have if we assumed that on the surface 
of the fluid there lay a solid crust capable of resisting 
tangential deformation, but opposing a negligible resistance 
to vertical deformation. 
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Application of the Fourier-Bessel inversion 
formula yields a first order differential equation 
for A as a function of ¢. 


2ndA /dt+ pgA/A= Jo” (da/dt)Jo(Ar)rdr. 
The solution is* 


(2.7) 


1 t 
f f Jo(Ar)rdrdt, (2.81) 
2n 


where K(A) must be determined from the initial 
conditions. Integrating the second term by parts 
with respect to ¢ 


‘do pg 
f —eratl 2nd t= geratl/2ar — f 
Ot 2nr 


Hence 


1 
4. ik oJ (dr)rdr 
2n vo 


1 
2n 0 


pg y 
f f Jo(dr)rdrdt. (2.82) 


This expression, with Eqs. (2.51), (2.52), (2.53) 
is the formal solution of the problem. 


3. SUBSIDENCE OF A CYLINDRICAL Bopy 


If we suppose that the fluid is initially at rest 
and that at ¢=0 a uniform circular load of radius 
ro is placed on the surface, 


0; r>Po 
o(r,t)=40; (3.1) 
o=const ; 
Then K(A)=0, and by using 
So"? Jo(dr)rdr= (ro/)Ji(Ar) 
A = ™ TJ, (Aro). (3.2) 
By substituting in (2.51), (2.52) and (2.53), 
(3.31) 
‘E. L. Ince, Ordinary Differential Equations, p. 21. 


r<ro, t>0. 
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V.=(oro/2n) 

XJIi(Aro) Jo(Ar)((1+Az)/A)dA, (3.32) 

p=oro fy J, (Aro) Jo(Ar)dr, (3.33) 
c= So'V.(r, 0, t)dt= (oro/ pg) fo°(1 — 

(3.34) 

We may note as a check that as [> 

pg) So” Ji(ro) Jo(Ar) dd 
0; 


(3.35) 


so that the system approaches the configuration 
of hydrostatic equilibrium asymptotically as it 
should. 

It does not appear to be possible to express the 
above integrals in closed form. They may be put 
in dimensionless form for numerical or mechan- 
ical integration by setting 


t=pgrot/2n, y=2/ro, 


t=p/o, pgi/o. 

Then 

vr=yfo (3.41) 


7/K-KuJ,(K) 

Jo(Kx)((14+Ky)/K)dK, (3.42) 
r= (3.43) 
w= (3.44) 


4. RECOVERY AFTER REMOVAL OF LOAD 


Setting the applied pressure, ¢, equal to zero 
over the whole surface in (2.82) we have 


A = 


K(X) can then be determined either from the 
initial velocity at the surface or from the initial 
configuration of the surface. From (2.52) 


V(r, 0,0)= (4.11) 


Hence, by inversion, 


K(A)=Af0rV .(r, 0, 0)Jo(Ar)rdr. (4.12) 


We also have 
0) + 0, 
=¢(r, 0) +(2n/pg) (A) (1 
XJo(Ar)add. (4.21) 


As ¢ becomes infinite, ¢ must approach zero, 
therefore 


0) = —(2n/pg) (4.22) 


or, inverting, 
= —(pg/2n) (r, 0)Jo(ar)rdr. (4.23) 


Thus the subsequent motion is completely deter- 
mined if we know either ¢(r, 0) or V,(r, 0, 0). 


5. APPLICATION TO Post-GLACIAL UPLIFT OF 
FENNOSCANDIA 


It would, of course, be impossible to calculate 
the motion of the earth throughout the history of 
the retreat of a continental ice sheet, since both 
the radius and the thickness varied with time in 
an irregular and, in the case of the thickness, an 
unknown manner. Moreover, an important part 
of the uplift (estimated as one-sixth to one-quar- 
ter of the total) was due not to plastic flow but 
to elastic recovery. However, ‘the results of the 
last section show that we may treat the part of 
the uplift which took place after the complete 
disappearance of the ice without any knowledge 
of how fast it retreated or how its thickness 
varied provided we know the rate of uplift at all 
distances from the center, or the configuration of 
the surface, at some time subsequent to the 
completion of the melting. During this part of the 
motion we do not need to consider the elastic 
component, since that would keep pace with the 
decrease in the load and would be complete when 
the ice had disappeared. 

Nansen’s curves (Figs. 1 and 2) give the 
amount and rate of uplift at a point in Angerman- 
land near the center of the Fennoscandian glaci- 
ated region and at Oslo, about 500 km distant. 
In order to use this data it will be necessary to 
assume a reasonable form for either ¢(r, 0), the 
initial configuration of the surface, or V2(r, 0, 0), 
the initial rate of uplift, and fit the assumed func- 
tion to the values given for these two points. 
It is reasonable to suppose that V,(r, 0, 0) will 
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Fic. 1. Uplift of the Fennoscandian region. Abscissas are 
times in units of 1000 years B.C. and A.D. Ordinates 
(solid curves) are depths of the surface in meters below 
equilibrium level. Ordinates of the dotted curve give the 
radius of the ice cap in kilometers. 


have a maximum at the center and will decrease 
outward in a way that can be represented with 
sufficient accuracy by an exponential function. 
We shall suppose then that 


V.(r, 0, 0) = —ae-”, (5.1) 


using the square of r in the exponent rather than 
the first power in order to avoid having a cusp 
at the center. By substituting in (4.12), 


K(A) = —an fore" Jo(Ar)rdr 
= 
By using (4.21) and (4.22), 
= (an/pgb*) fore! Ndr. (5.3) 
At r=0, t=0 
£(0, 0) = (an/pgb?) = (nab/ pg). 
The kinematic viscosity, y= 7/p, is then given by 
v=gf(0, 0)/2./rab. (5.41) 


Since we are at liberty, in applying the above 
theory, to take any time after the disappearance 
of the ice as the initial instant ‘= 0, we may check 
roughly the validity of the assumption that 
V.(r, 0, 0) does not depart very widely from the 
form given in (5.1) by calculating v from (5.41) 
for various initial instants and seeing that it re- 
mains of the same order of magnitude. If we ex- 
press ¢(0, 0) in meters, a in meters per century, 
and 6 in reciprocal kilometers, v is given in c.g.s. 
units by 


(5.2) 


v=8.73 X10"¢(0, 0)/ad. (5.42) 
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Fic. 2. Rate of uplift in meters per century. Abscissas in 
units of 1000 years. 


¢(0, 0) may then be read from the curve for 
Angermanland in Fig. 1, a from that in Fig. 2, 
and 5 may be calculated from the curve for Oslo 
in Fig. 2 by setting r= 500 in (5.1). In Table I the 
values of these quantities are given for various 
times taken as the initial instant. 

The initial configuration of the surface, corre- 
sponding to the assumed initial velocity (5.1), 
is given by 


¢(r, 0) = (an/pgb?) fore” Jo(Ar) (5.51) 
or 
0)/¢(0, 0) 

= (1/21) 4 (5.52) 


where x= br. Values of this quantity as a function 
of x are given in Table II. Taking the initial 
instant at 5000 B.C., and using b=1.27 x10" 
km", the initial profile of the surface and rate of 
uplift are plotted in Fig. 3. A complete observed 
profile is not available to compare with the calcu- 
lated curve in order to determine how closely the 
assumptions made fit the actual case, but it may 
be noted that Fig. 3 indicates a depression of the 
surface at 500 km radius of about 76 meters. 
This is to be compared with 85 meters at Oslo at 
5000 B.C. as read from the curve of Fig. 1. 

In comparing the figure given above for the 
earth’s kinematic viscosity with experimentally 
determined viscosities of solid substances it 
should be noted that two quite distinct properties 


TABLE I. Values of the quantities entering in Eq. (5.42) for 
various times taken as the initial instant. 


t=0 a bx10® 
5,000 B.C 3.9 1.27+.07 147 2.6+0.2 
4,000 “ 2.7 1.19+.09 118 3.240.3 
3,000 “ y 1.24+.11 94 3.0+0.4 
2,000 « ‘ 1.8 1.25+.15 74 2.940.5 
Mean value: »=2.9X102! cm? sec.~ 
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(50,0 Vz (0,0,0) = -3.9 mpececentvry 
Vz, (0,0 0) 
2000 1509 1000 500 1000 1700 2000 
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Fic. 3. Rate of uplift and surface profile at 5000 B.C. Vertical exaggeration of profile 1.7 <10* times. 


TABLE II. 

x 0)/¢(0, 0) x 0)/¢(0, 0) 
0.0 +1.000 1.0 +0.155 
0.2 +0.922 1.5 —0.076 
0.4 +0.780 2.0 . —0.068 
0.6 +0.562 —0.033 
0.8 +0.328 


of matter are often included in the term.® One, 
which is better designated as inner friction, relates 
to the damping of elastic vibrations, the other 
relates to plastic flow. The coefficient of inner 
friction is of the order of 10° or 10° c.g.s. units for 
most metals and is higher for soft than for hard 
substances. The ‘‘true”’ coefficient of viscosity is 
more difficult to measure experimentally, differ- 
ent observers giving widely varying values,® and, 
in general, increases rapidly with the hardness of 
the substance. The figures given in Table III are 
taken from Gutenberg and Schlechtweg’s paper. 
In view of the great increase of viscosity with the 
hardness of the material and the pressure, it is not 
surprising to find that for silicates under the 
enormous pressures within the earth the viscosity 
is as high as we have found here, even when the 
decrease of viscosity at high temperatures is 
taken into consideration. 

R. W. van Bemmelen and H. P. Berlage’ have 
recently made an estimate of the viscosity of the 
earth on the basis of post-glacial uplift using a 


5B. Gutenberg and H. S. Schlechtweg, Viscosity and 
Inner Friction of Solid Bodies, Physik. Zeits. 31, 745 (1930). 

*W. D. Kusnezow, Zeits. f. Physik. 51, 239 (1928). 

7R. W. van Bemmelen and H. P. Berlage, Gerl. Beit. z. 
Geophys. 43, 19 (1934). 


TABLE III. Values taken for the coefficient of viscosity of 
various substances. 


MATERIAL TEMPERATURE (g. cm™~! sec.~!) 
Shoemakers wax 2 108 
Asphalt 15 3 
Ice 0 5 x10" 
Glass 575 1 x10" 
Lead 20 1 x10" 
Calcite 18 1.510" 
Rocksalt 18 2 


method quite different from that followed hree. 
They find »= 1.3 X10” which, since they assume 
a density of 3, is equivalent to y= 0.4 X 10°. Their 
method makes use of special hypotheses concern- 
ing the structure and dynamics of the earth’s 
crust and is only indirectly based on hydro- 
dynamical theory. They also assume that the 
flow is confined to a layer 100 km thick and it is 
therefore to be expected that their value should 
be lower than that found by the present method. 
According to their theory v varies as the cube of 
the assumed thickness of the layer, hence if it 
were taken to be 400 km, their figure would be 
comparable with ours. 

In conclusion the author wishes to express his 
thanks to Professor L. B. Slichter, under whose 
guidance this work was done, for many helpful 
suggestions and criticisms. The graphical evalua- 
tion of the integral given in Table II was very 
kindly carried out by members of the Mathe- 
matics Laboratory of the Massachusetts Institute 
of Technology. Acknowledgment should also be 
made to Professor R. A. Daly’s stimulating lec- 
tures on Pleistocene movements of the earth’s 
crust, which originally suggested the problem. 
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A Note on Viscosity as a Function of Volume and Temperature of Oils 


. R. B. Dow, Research Laboratory of Physics, Harvard University 
(Received April 17, 1935) 


The viscosity-volume data of Kleinschmidt and Dow 
have been examined at various pressures and temperatures 
for lard, sperm and Pennsylvania medium oil, The viscosity- 
volume isotherms at 25°, 40° and 75°C are not identical for 
any of the oils studied, indicating that viscosity cannot be 
a function of the specific volume alone. The viscosity- 
volume curve for lard oil at 25° departs from the one at 75° 
by an amount sufficient to change the viscosity by a 


factor of 2.3 at a volume of 0.99, and by a factor of 3.2 ata 
volume of 0.93. Similar curves for Pennsylvania medium 
oil at the same temperatures are even more relatively dis- 
placed; the discrepancy in viscosity varies from a factor 
of 3.8 at a volume of 0.99 to 7.6 at 0.94. The three oils do 
not obey Batschinski’s equation at atmospheric and higher 
pressures up to 4000 kg/cm?*. 


OMPARATIVELY little is known of the 
physical properties of lubricating oils at 
high hydrostatic pressures. Among the non- 
thermodynamic properties of lubricants at high 
pressures, viscosity has been most extensively 
studied because of its significance for thick film 
lubrication. The experiments of Hyde! and, more 
recently, those of Hersey and Shore,” and Klein- 
schmidt,? have shown that the coefficient of 
viscosity of a mineral oil at ordinary temperatures 
increases by a factor of about 20 with an initial 
increase of pressure of 1000 kg/cm,” this increase 
being several times greater than that observed 
for pure liquids‘ or mixtures of liquids’ through 
the same range of pressure. With the recent 
study of some of the thermodynamic properties 
of similar oils,* sufficient data are available for an 
examination of the viscosity of oils as a function 
of volume. 

In addition to the practical usefulness of 
viscosity-volume data taken at various pressures 
and temperatures, there is theoretical interest in 
the functional relation between viscosity and 
volume. Consequently, this communication pre- 
sents the viscosity-volume-temperature relations 
for three lubricating oils and includes a discussion 
of the theoretical relationship. 


DATA 


Table I contains the log relative viscosities and 
volumes at various pressures and temperatures 


a a Hyde, Proc. Roy. Soc. A97, 240 (1920). 

2M. somal and H. Shore, Mech. Eng. 50, 221 (1928). 
3 R. V. Kleinschmidt, Trans. A.S.M.E. APM-50-4 (1928). 
4P. W. Bridgman, Proc. Am. Acad. 61, 57 (1926). 

5 R. B. Dow, Physics 6, 71 (1935). 

*R. B. Dow, J. Wash. Acad. Sci. 24, 516 (1934). 


for lard, sperm and Pennsylvania medium oil, 
respectively, the data being taken from the 
papers of Kleinschmidt* and Dow.* The density 
of each oil at atmospheric pressure and 40°C js 
given in order that the specific volumes may be 
computed directly from the table of volumes by 
division. The log relative viscosities are expressed 
as logio t/to, ¢ being the time of fall of a weight in 
a viscometer at a certain pressure and tempera- 


TABLE I. Relative viscosity and volume. 


logio RELATIVE VISCOSITY VOLUME 
PRESSURE 
KG/cM? 25° 40° 75° 25° 40° 75° 
Lard oil pw =0.9009 g/cm? 
1 0 1.770 1.370 0.9902 1.0000 1.0190 
100 0.079 .845 -9936 
250 -183 .938 500 .9763 -9850 1.0051 
500 345 0.082 628 -9647 -9721 -9921 
750 499 .220 -742 -9550 .9800 
1000 .642 351 -9523 -9697 
1500 .920 .607 0.070 9299 9366 -9522 
2000 .835 .262 9229 
2500 1.052 441 9111 9240 
3000 615 -9120 
4000 962 .8927 
Sperm oil pw =0.8945 
1 0 1.720 1.256 0.9894 1.0000 1.0227 
100 0.079 -792 -9835 
200 -150 1.0099 
300 .220 9818 
400 .289 981 -9685 9768 
500 0.040 -9722 
750 .9618 .9794 
1000 .318 .757 .9684 
1500 959 .9437 -9510 
2000 0.149 -9368 
2500 .327 9241 
3000 481 -9127 
4000 792 .8926 
Pennsylvania oil po =0.8524 g/cm® 
1 0 1.660 1.020* 0.9901 1.0000 1.0178 
100 0.119 .761 .9839 
250 .280 .235* -9752 9841 1.0040 
500 536 0.131 .420* -9632 9711 .9908 
750 777 .594* .9599 .9786 
1000 1.008 551 -760* -.9440 -9504 .9672 
1500 955 0.070 9340 .9485 
2000 1.341 .369 9196 .9330 
2500 
3000 -953 -9082 
4000 1.511 


* Extrapolated. 
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Fic. 1. Relative viscosity of lard oil as a function of volume. 


ture and ¢) the time of fall at atmospheric pres- 
sure at 25°. The reader is referred to the original 
papers for details of method and computation. 


DISCUSSION 


It is apparent at once on inspection of Figs. 1, 
2, and 3, which are drawn from the data of 
Table I, that the relative viscosity increases 
greatly with comparatively small decrease of 
volume. Inasmuch as the figures are drawn from 
data of different observers, as mentioned previ- 
ously, there is some doubt as to the experimental 
accuracy of the viscosity-volume curves. How- 
ever, the writer’s acquaintance with the method 
used by Kleinschmidt for the viscosity deter- 
minations leads him to estimate the inaccuracy 
of the curves to be not more than a few percent, 
which does not seriously limit their applicability. 
The lack of serious deviation of the points repre- 
senting the experimental values from the smooth 
curves of the figures gives an indication of the 
probable degree of accuracy of the data. 

An interesting feature of the figures is the 
relative displacements of the viscosity-volume 
curves at 25°, 40° and 75°C. If viscosity were a 
function of volume only, the curves for each oil 
would coincide at all three temperatures. The 
figures show that this requirement is not satisfied 
in any case. The viscosity-volume curve for lard 
oil at 25° departs from the curve at 75° by an 
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Fic. 2. Relative viscosity of sperm oil as a function of 
volume. 
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Fic. 3. Relative viscosity of Pennsylvania medium oil as 
a function of volume. 


amount sufficient to change the viscosity by a 
factor of 2.3 at a volume of 0.99, and by a factor 
of 3.2 at a volume of 0.93. Similar curves for the 
Pennsylvania oil at 25° and 75° are even more 
relatively displaced; the discrepancy in viscosity 
varies from a factor of 3.8 at a volume of 0.99 to 
7.6 at 0.94. 


R. B. 


Of the theoretical relationships that have been 
proposed for viscosity as a function of volume, 
that of Batschinski? has been the most useful. 
His equation states that viscosity is related to 
the specific volume in the following way 


n=c/v—w. 


c and w are characteristic constants of the liquid. 
Batschinski showed that if the change of fluidity 
of normal liquids, such as benzene and ethy! 
ether, was expressed with change of volume 
caused by change of either temperature or pres- 
‘sure, a linear relation resulted as demanded by 
his equation. But the data on viscosity and volume 
as functions of pressure were so limited in range of 
pressure available at that time, that both quanti- 
ties could be fairly well expressed as linear func- 
tions of pressure. It is well known that such func- 
tions do not remain linear through the range of 
pressure now available. 

Batschinski’s equation has been used lately by 
Bingham and Brown,’ Bingham and Coombs,’ 
and Lederer! in deducing theories of viscosity. 
Although there is little experimental information 
concerning lubricants in this respect, R. N. J. 
Saal" in discussing the influence of pressure on the 
viscosity of a nonplastic, asphaltic bitumen, con- 
sidered that his experimental data showed that 
the decrease of viscosity as the temperature rose 
was due to the effect of thermal expansion. Thus 
in addition to what has already been established 
without any assumption as to theory, namely, 
that viscosity is not a function of volume only, 
it is desirable to call attention to the failure of 
Batschinski’ equation at high pressures, for this 
limitation does not seem to have been generaly 
established. 

Bridgman* has published viscosity-volume 
data for several normal liquids. His results gave 


7 A. Batschinski, Zeits. f. physik. Chemie 84, 643 (1913). 
5’ E. C. Bingham and D. F. Brown, J. Rheol. 3, 95 (1932). 
®E. C. Bingham and C. E. Coombs, Physics (New York 
Meeting). 
E. L. Lederer, Kolloid-Bei. 34-35, 270 (1932). 
(1934) N. J. Saal, Proc. Wor. Pet. Congress, London, 521 
12 P, W. Bridgman, Proc. Am. Acad. 66, 185 (1931). 
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curves similar to the figures of this paper. An 
examination of Batschinski’s equation by substi- 
tution of Bridgman’s data gives a general result 
that in no case is a linear relation obtained be- 
tween fluidity and specific volume. Benzene, 
ethyl ether, pentane, etc., obey the Batschinskj 
equation at atmospheric pressure, that is when 
the equation is applied for change of temperature, 
but at higher pressures the invalidity of the equa- 
tion is beyond experimental error. 

Consequently, on referring to Figs. 1, 2, and 3 
again, it is not unexpected that the fluidity curves 
would not bear a linear relation to the volume, 
and the observed displacements of the curves 
show that the viscosity, or fluidity, is also a func- 
tion of temperature. Accordingly. the constants ¢ 
and w of Batschinski’s equation vary with pres- 
sure and temperature for these three oils. These 
data establish a point of much theoretical inter- 
est, namely that pressure and temperature 
changes affect viscosity differently. 

Bridgman" in a recent paper on some of the 
theoretical aspects of high pressure phenomena 
has discussed the effects of temperature and 
pressure on the energy of solids, showing that in 
the case of NaCl the change of energy internal to 
the atom is nearly three times as great when a 
definite change of volume is brought about by a 
change of pressure as when brought about by a 
change of temperature. His explanation considers 
a compressible atom as demanded by a theorem 
of Schottky. It is likely that the fundamentals of 
the situation apply to liquids and give a possible 
explanation of the different effects of pressure 
and temperature on viscosity, although the 
problem remains to be treated quantitatively. 
The experimental evidence for the compressible 
atom is uncontrovertible to such an extent that 
Batschinski’s conception of atomic volume con- 
stants cannot be valid over a wide experimental 
range. It is the desire of the author to call atten- 
tion to these serious limitations of Batschinski’s 
theory as demanded by the experimental data, 
rather than to question. the usefulness of the rela- 
tion at atmospheric pressure. 


8 P, W. Bridgman, Rev. Mod. Phys. 7, 6 (1935). 
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Since in a colloidal or quasi-plastic material, the viscosity 
is a function of the rate of shear, it seems best to determine 
the velocity-gradient across a section perpendicular to the 
direction of flow under shearing forces, and to determine 
the coefficient of viscosity under a given velocity gradient. 
An account is given of the experimental methods at present 


available for this purpose, and a method indicated by 
which the viscosity may be calculated from such measure- 
ments. The viscometers described comprise instruments 
based on the optical, electromagnetic, hot-wire and color- 
band principles, respectively. 


HAT a colloid behaves abnormally under 
shearing stresses has been known since 
Bingham! first investigated the matter. As a 
consequence of this fact, any attempt to deduce 
a coefficient of viscosity for such a substance from 
measurements in either a capillary-flow or a 
concentric-cylinder apparatus will only produce 
qualitative results, for the viscosity is a function 
of the rate of shear, which will vary across a 
section of the flow. It was Hatschek? who sug- 
gested that if the velocity-gradient across a sec- 
tion of the fluid in motion could be measured, it 
would be possible to determine the coefficient of 
viscosity in terms of the velocity-gradient. In 
spite of ingenious attempts to assume some such 
function and to substitute it in the Poiseuille or 
Couette equations, as the case may be, and de- 
duce the form of this function from the usual 
measurements of outflow or torque, respectively, 
it is not until one actually measures the velocity 
or its gradient from point to point across a section 
of the flow that one can calculate the viscosity as 
a function of the shear with certainty. The object 
of the present paper is to describe the methods at 
present available. They are all of quite recent 
origin. They may be grouped under four heads, 
viz. : 
(1) The photographic method. Particles moving 
in the fluid reveal their relative velocity by the 
relative lengths of the streaks which their images 
make on a photographic plate during a certain 
exposure. 

(2) The electromagnetic method. This depends 
on the induction produced by charged particles 
in the fluid moving in a magnetic field. 

(3) The hot-wire method. The cooling of the wire 


' Bingham, Fluidity and Plasticity (1922), New York. 
* Hatschek, Viscosity of Liquids (1928), London. 


is proportioned to the velocity at that point in 
the fluid. 

(4) The visual method. By a color scheme im- 
posed on the substance, it is made to exhibit the 
form into which an originally plane section is 
curved under shearing stress. 

The photographic method was apparently first 
used by Porter and Nisi* in 1923 for studying the 
flow of air round obstacles on a small scale. They 
placed a ball bearing 3 mm (1/8’’) diameter in a 
tube of square section at a point where a thin 
beam of light could be directed across the tube. 
To make the motion visible smoke was let into the 
tube at the level of this illuminated plane, and 
the camera directed in a direction perpendicular 
to it. When an exposure is made, the particles of 
smoke describe streaks on the photographic plate, 
whose lengths are proportional to their relative 
velocities. The method in some modification or 
other has been frequently employed in hydro- 
dynamics since that time, though not often for 
the study of colloidal flow. In a sense, of course, 
the medium employed by Porter and his collabo- 
rators was colloidal, but the admixture of smoke 
was not at sufficiently high concentration to 
produce abnormal flow. The success of the 
method lies in the obtaining of the thin pencil of 
very brilliant light and in using the scattered 
light from the solid phase rather than the light 
directly transmitted through the tube. Pichot 
and Dupin‘ have used this method on a fluid par- 
taking of both suspension and sol characteristics, 
consisting as it did of a7 percent solution of gela- 
tine in water with aluminum particles added to 
make visible the relative velocities of different 
parts of the fluid in a tube when an exposure was 


3 Porter and Nisi, Phil. Mag. 46, 754 (1923). 
4 Pichot and Dupin, Comptes rendus 192, 1079 (1931). 
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made. From the only two photographs published 
by them (showing water and gelatine solution 
severally in the same tube) it is difficult to detect 
any considerable difference in the behavior of 
homogeneous and heterogeneous liquids. 

By Faraday’s laws of electromagnetic induc- 
tion, when a conducting liquid is in motion in a 
magnetic field, an electromotive force is induced 
in it whose value is proportional to the product of 
the magnetic force and the velocity. Williams’ at- 
tempted to measure the velocity-gradient across 
a tube in which an electrolyte was flowing by ap- 
plying a magnetic field along a diameter and 
measuring the gradient of electromotive force 
along a perpendicular diameter on a sensitive 
galvanometer. Actually the instrument measures 
the voltage between two points and it is necessary 
to take one of these points on the interior surface 
of the tube, while the other is a probe, insulated 
save for the point, which is traversed along a 
diameter. The experimental tube was 1 cm diam- 
eter and the liquid, copper sulphate solution. A 
comparison is made of the theoretical voltage due 
to Poiseuille flow for a given distance of separa- 
tion of the electrodes (expressed as a fraction of 
the diameter) and the observed voltage. It is not 
clear, however, how one could readily work back 
from a series of measured voltages to the un- 
known velocity distribution. With a string 
galvanometer, the method would serve for 
rapidly fluctuating velocities as in turbulent 
motion. 

It is of interest to note that the method was 
used earlier on a large scale investigation of the 
tidal motion in the estuary of the River Dart 
(Devon). 

We may conveniently include in this section 
the method of Alcock and Sadron,‘ since it is 
based on Maxwell’s electromagnetic theory of 
light and was demonstrated by him in Canada 
balsam. A narrow beam of polarized light is 
doubly refracted by its passage through a gypsum 
plate. This refraction is modified by the subse- 
quent transmission of the light through a stratum 
of the moving fluid, causing a shift of the inter- 
ference bands arising from the double refraction 
to an extent proportional to the velocity of the 
stratum. As this method has been recently de- 


5 Williams, Proc. Phys. Soc. London 42, 466 (1930). 
6 Alcock and Sadron, Physics 6, 92 (1935). 
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Fic. 1. Hot-wire anemometer. 


scribed in detail in this Journal, it is not proposed 
to repeat the description now. 

An instrument which is still widely used in 
aerodynamics for measuring velocity is the Pitot 
tube. This consists essentially of a tube having a 
fine orifice facing upstream. The pressure differ- 
ence between the air in this tube and in another 
shielded from the direct current (the static tube) 
is recorded on a sensitive manometer and is by 
Torricelli’s law proportional to the square of the 
velocity. Kroepelin’ has traversed such an ane- 
mometer across the end of a tube in which a 
gelatine solution was flowing at the place where 
it debouched upon the outflow reservoir and re- 
corded small departures from the Poiseuille 
(parabolic) distribution. 

A more sensitive anemometer for these pur- 
poses is that which depends upon the cooling 
suffered by a wire, electrically heated to about 
25 degrees above the fluid temperature, when the 
fluid flows past it. The cooling is measured on a 
Wheatstone bridge in terms of the change of re- 
sistance experienced by the wire. Before use, the 


7 Kroepelin, Zeits. f. physik. Chemie A149, 291 (1930). 
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Fic. 2. Typical velocity-gradients in a suspension of 
starch in a solvent mixture of carbon tetrachloride and 


paraffin for four different rates of rotation of the outer 
cylinder of the hot-wire anemometer. 


hot-wire anemometer must be calibrated by 
carrying it on a small whirling arm through the 
same still fluid contained in a tank. The calibra- 
tion curve is obtained in the form of resistance as 
a function of the velocity. It is then ready to be 
applied to the measurement of velocity-gradients 
in the fluid in question. Instruments in which 
liquid flows through an oval tube, the wire being 
parallel to the flats of the oval section, or of the 
concentric cylinder type in which the wire is 
vertical and parallel to the axis of the cylinders, 
have been constructed by the writer, and one of 
the latter (Couette) type has been used by Dr. 
Tyler® and the writer? with some success for the 
study of velocity-gradients in suspensions. 

The apparatus is shown in Fig. 1. The hot-wire 
W is of nickel 0.025 mm (0.001’’) diameter and 
1.5 cm (0.6’’) long and is heated by a current of 
0.3 amp. It is soldered to the ends of two fine 
needles inserted through ebonite sleeves in the 


* Richardson and Tyler, Proc. Phys. Soc. London 44, 
142 (1933). 


* Richardson, Les A ppareils a Fil Chaud (1934), Paris. 


METHODS 275 
outer cylinder which is rotated through the pulley 
P below. Electrical connection to the needles is 
assured by dippers working in the fixed annular 
rings of mercury GG mounted below the fork. 
The speed of this cylinder is observed by means 
of a Neon lamp N giving 50 flashes per sec. which 
illuminates the stroboscopic disk D. The hot- 
wire can be traversed across the annular space 
between the cylinders which are, respectively, 
1.27 and 2.54 cm in diameter by means of the 
micrometer head M. The torque communicated 
to the inner cylinder may be taken up by a phos- 
phor bronze suspension wire (not shown in the 
sketch) the twist on which was indicated by the 
usual mirror, lamp and scale. Fig. 2 shows typical 
velocity-gradients in a suspension of starch in a 
solvent mixture of carbon tetrachloride and paraf- 
fin (54.5 percent) which had the consistency of a 
paste (circles or crosses and continuous lines) for 
four different rates of rotation of the outer 
cylinder. On the same figure is shown one set 
(dots and broken line) for the solvent alone. No- 
tice that although the latter, in accordance with 
theory, exhibits a nearly constant velocity- 
gradient through the annular space, the former 
shows a gradient growing steeper towards the 
inner cylinder, which is located on the left-hand 
ordinate of the figure. These curves are typical of 
the observations in starch suspensions as well as 
of clay pastes and gelatine solutions. 

By using an Einthoven string galvanometer 
coupled through a transformer to the hot-wire 
circuit, it is possible to record the onset of turbu- 
lence in fluid flow at any point by noticing when 
the current through the wire begins to oscillate 
due to fluctuations in the heat losses from it 
engendered by turbulence. 

Coming finally to the method by which the 
velocity-gradient is made visible by a dye, we 
find in a recent paper by Lawrence"® application 
of this device to ammonium oleate solutions, 
known to show anomalous flow. By means of a 
special two-way cock, colored and uncolored fluid 
were alternately admitted to a tube. The dye used 
was the ammonium salt of Congo red. The head of 
the dye on entering the tube had a plane surface 
but after it had proceeded a little way into the 
flow tube it was distorted into a form which de- 
lineated the velocity-gradient. An exposure was 


10 Lawrence, Proc. Roy. Soc. A148, 59 (1935). 
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Fic. 3. Photograph of one of Pearson’s metal billets extruded through a die. 


then made and the resulting photograph, when 
enlarged, enabled the required velocity-gradient 
to be measured up. In a homogeneous fluid the 
heads of the colored sections would, of course, 
have been parabolas, but in the colloids devia- 
tions from these curves were obtained. 

As another instance of the application of the 
visual method may be cited the use of a painted 
or etched pattern in the study of the extrusion of 
ceramic materials or metals through dies. This 
method seems to have been first applied by Trit- 
ton" to plasticine, alternate layers of two colors 
of which were put in a press and forced out 
through a hole. Mr. Pearson of the Metallurgical 
Department of this University has studied the 
extrusion of alloys by such a method, although 
being solid much greater pressures are required 
than with any of the other media we have men- 
tioned. Pearson took a “billet” of the metal, cut 
it down the axis, painted the exposed plane sur- 
face of one of the half cylinders thus formed with 
ink in transverse stripes and put the two halves 
together again. The reunited billet was then 
placed in the barrel of a press which it just fitted, 
and the metal extruded through a die of aperture 
1 cm diameter at the bottom of the press. As the 
alloy “flowed out,” the stripes were distorted into 
shapes which represented the velocity-gradients 
under the given ‘‘head”’ of pressure. The pattern 
was exposed by separating again the two halves 
of the billet. A photograph of one of Pearson’s 
extruded billets is shown (on its side) in Fig. 3. 

In attempting a critique of the various methods 
that he has described above, the author is faced 
with the difficulty that with only two methods— 
and one of them his own—is a systematic meas- 
urement of a series of gradients available. The 
remaining workers have been content to show 


" Tritton, J. Inst. Met. 26, 250 (1921). 


isolated specimens of the capabilities of their de- 
vices. There are, however, certain directions in 
which comparisons may be useful. One question 
which the colloid physicist will be sure to ask is 
as to the extent with which the measuring devices 
will interfere with the flow to be observed. Thus, 
in applying the photographic method one may 
have to introduce strongly light-scattering par- 
ticles into the medium, which may change its 
properties. On the other hand this same method 
is well adapted to the study of systems which by 
their very nature contain solid particles, such as 
paints. In the electrical and hot-wire methods the 
electrodes may cause some interference, but since 
they may be formed of quite fine wires this source 
of trouble is not prohibitive. To the extent that 
the hot-wire is above the temperature of the 


surrounding fluid it may cause local heating and’ 


convection currents, but the excess of tempera- 
ture is not more than 25 degrees and need not 
exceed 10 degrees. The chief difficulty in the 
color method appears to lie in the discontinuity 
in the flow at the change-over from one supply to 
the other, although this can be obviated by a 
well-designed two-way cock. The method of 
Lawrence requires the original colloid to be trans- 
parent, as does also that of Alcock and Sadron, 
while the method of Williams is restricted to 
electrolytes. Many of the authors have tested 
their respective instruments on an homogeneous 
fluid in streamline flow, and this is certainly a 
good test of the freedom of the apparatus from 
interference to the flow, which would at once be 
apparent as a deviation from parabolic distri- 
bution. 

As typical of the types of flow recorded in tubes 
below the critical velocity we have grouped three 
curves on Fig. 4: (A) Parabolic distribution in 
paraffin, recorded by hot-wire, (B) a somewhat 
flattened curve in ammonium oleate, recorded by 
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Fic. 4. Velocity distributions in tube flow. A, parabolic 


_distribution in paraffin; B, a somewhat flattened curve in 


ammonium oleate; C, ‘‘plug flow” in extruded plasticine. 


Lawrence’s dye method, (C) “plug flow’’ in 
extruded plasticine, shown by the shape of the 
head of the emergent stem. 

It will probably be asked how these data may 
be used to measure the viscosity of the substance. 
To do this it is necessary to obtain a formula for 
the velocity-gradient and exhibit how the vis- 
cosity varies with this by means of a graph. Thus 
if the apparatus is of the flow-tube type, the 
forces in equilibrium on a cylinder of radius r and 
of unit length are 


ar p=2nrn-dv/dr, 


where p is the pressure inducing flow and dv/dr 
is the velocity gradient and 7 the coefficient of 
viscosity. Therefore n is proportional to pr/dv/dr 
and we only need to plot this latter quantity 
against dv/dr to exhibit the way in which the 
viscosity varies with gradient. 

In the case of the concentric cylinders, the 
velocity gradient at any radial distance being 
rdw/dr, where w is the angular velocity of the 
outer cylinder, the moment of the annulus about 
the axis is 2rr*ndw/dr per unit length, and this 
moment is constant for all annuli, so that the 
torque is proportional to r*ndw/dr. Thus 
t/(r'dw/dr) should be a measure of the relative 
viscosity. Thus from Fig. 2, the slope of the 
angular velocity curves may be plotted against r, 
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Fic. 5. Variation of viscosity with velocity-gradient. 


and from this another curve of r/(r’dw/dr) could 
be plotted against rdw/dr, r being known from the 
twist on the inner cylinder. This curve is shown 
on Fig. 5 and accordingly shows us how the vis- 
cosity varies with velocity-gradient. The hori- 
zontal line shows the viscosity of the solvent. 

Such methods as described in this paper may 
become more common in rheology in the near 
future, and may lead to the coefficient of viscosity 
for colloidal substances being defined at a 
standard velocity-gradient, for example 1 cm/sec. 
per cm. This is certainly preferable to the meth- 
ods by which substances are put successively into 
a viscometer without any attempt to see that 
they are tested at comparable velocity-gradients, 
so that the two substances are being compared 
at varying and quite unknown rates of shear. 
Although technical viscometers are favored by 
industrial chemists on account of their simplicity 
and rapidity of working, it is to be feared that 
any results they may obtain with a substance 
whose viscosity varies in a fashion similar to that 
shown in Fig. 5—and this class includes a great 
many industrial products—must be quite value- 
less. 


